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Abstract 



The goal of this paper is to study the lookdown model with selection in the case 
of a population containing two types of individuals, with a reproduction model which 
is dual to the A-coalescent. In particular we formulate the infinite population "A- 
lookdown model with selection". When the measure A gives no mass to 0, we show 
that the proportion of one of the two types converges, as the population size N tends to 
infinity, towards the solution of a stochastic differential equation driven by a Poisson 
point process. We show that one of the two types fixates in finite time if and only 
if the A-coalescent comes down from infinity. We also consider the general case of a 
combination of the Kingman and the A-lookdown model. 

Subject classification 60G09, 60H10, 92D25. 

Keywords Look-down with selection. Lambda coalescent, Fixation and non fixa- 
tion. 

1 Introduction 

In this paper we consider the lookdown (which is in fact usually called the "modified 
lookdown") model with selection when we replace the usual reproduction model by a 
population model dual to the A-coalescent. We first recall the models from [15] and [9], 
and then we will describe the variant which will be the subject of the present paper. 

Pitman |15j and Sagitov |16j have pointed at an important class of exchangeable 
coalescents whose laws can be characterized by an arbitrary finite measure A on [0, 1]. 
Specifically, a A-coalescent is a Markov process (Ilt,t > 0) on Voo (the set of partition 
of N) started from the partition Oqo := {1)2,...} and such that, for each integer 
n > 2, its restriction (n["', t > 0) to Vn (the set of the partition of {1, 2, ... , n}) is a 
continuous time Markov chain that evolves by coalescence events, and whose evolution 
can be described as follows. 
Consider the rates 




2<i<k. 



(1.1) 
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Starting from a partition in Vn with k non-empty blocks, for each £ = 2, . . . , /c, every 
possible merging of I blocks (the other k — t blocks remaining unchanged) occurs at rate 
\k,i, and no other transition is possible. This description of the restricted processes 
n" determines the law of the A-coalescent IT. 

Note that if A({0}) = A([0, 1]) > 0, then only pairwise merging occur, and the cor- 
responding A-coalescent is just a time rescaling (by A(0)) of the Kingman coalescent. 
When A({0}) = which we will assume except in the very last section of this paper, 
a realization of the A-coalescent can be constructed (as in [15]) using a Poisson point 
process 

oo 
i=l 

on M-i- X (0, 1] with intensity measure dt®v{dp) where u{dp) = p^'^A{dp). The measure 
v{dp) may have infinite total mass. Each atom (t,p) oi p influences the evolution as 
follows : 

• for each block of n(t^) run an independent Bernoulli {p) random variable; 

• all the blocks for which the Bernoulli outcome equals to 1 merge immediately 
into one single block, while all the other blocks remain unchanged. 

In order to obtain a construction for a general measure A, one can superimpose onto 
the A-coalescent independent pairwise merges at rate A({0}). 

The lookdown construction was first introduced by Donnely and Kurtz in 1996 [9]. 
Their goal was to give a construction of the Fleming- Viot superprocess that provides 
an explicit description of the genealogy of the individuals in a population. Donnelly 
and Kurtz subsequently modified their construction in [lOj to include more general 
measure-valued processes. Those authors extended their construction to the selective 
and recombination case |11] . 

We are going to present our model which we call A-lookdown model with selection. 
An important feature of our model is that we will describe it for a population of infinite 
size, thus retaining the great power of the lookdown cosntruction. As far as we know, 
this has not yet been done in the case of models with selection except in our previous 
publication [3], where we considered a model dual to Kingman's coalescent. 

We consider the case of two alleles b and B, where B has a selective advantage over 
b. This selective advantage is modelled by a death rate a for the type b individuals. 
We will consider the proportion of b individuals. The type b individuals are coded by 1, 
and the type B individuals by 0. We assume that the individuals are placed at time 
on levels 1,2,..., each one being, independently from the others, 1 with probability x, 
with probability 1 — x, for some < x < 1. For each i>l and t > 0, let rit{i) G {0, 1} 
denote the type of the individual sitting on level i at time t. The evolution of {r]t{i))i>i 
is governed by the two following mechanisms. 

1. Births Each atom {t,p) of the Poisson point process m corresponds to a birth 
event. To each {t,p) £ m, we associate a sequence of i.i.d Bernoulli random 
variables {Zi,i > 1) with parameter p. Let 

hp = {i>l: Zi = l}. 
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and 

lt,p = inf{z G It,p ■ i > min It^p} 

At time t, those levels with Zi=l and i > it^p modify their label to r/^-(min It,p)- 
In other words, each level in It^p immediately adopts the type of the smallest level 
participating in this birth event. For the remaining levels reassign the types so 
that their relative order immediately prior to this birth event is preserved. More 
precisely 



'rit-{i), if i < lt,p 

r]t- (min It^p), if i G It,p \ {min It^p} 

.Vt-ii - {#{It,p n [1, . . . ,z]} - 1)), otherwise 



We refer to the set It^p as a multi-arrow at time t, originating from min It^p, 
and with tips at all other points of It^p. This procedure is usually referred to 
as the modified lookdown construction of Donnelly and Kurtz. In the original 
construction, the types of the levels in the complement of It^p remained unchanged 
at time t, hence the types r)t-{i), for i G It,p \ {min It^p} got erased from the 
population at time t. 

2. Deaths Any type 1 individual dies at rate a, his vacant level being occupied by his 
right neighbor, who himself is replaced by his right neighbor, etc. In other words, 
independently of the above arrows, crosses are placed on all levels according to 
mutually independent rate a Poisson processes. Suppose there is a cross at level 
i at time t. If = 0, nothing happens. If r]f-{i) = 1, then 



Vt-{k), if k<i; 
rit-{k + l), if k>i. 



We refer the reader to Figure 1 for a pictural representation of our model. Note that 
the type of the newborn individuals are found by "looking down", while the type of 
the individual who replaces a dead is found by looking up. So maybe our model could 
be called "look-down, look-up". 

Since we have modelled selection by death events, the evolution of the N first 
individuals r]t{l) , . . . , rjtiN) depends upon the next ones, and X^^ = N^^{r]t{l) + 
...r]t{N)) is not a Markov process. We will show however that for each t > the 
collection of r.v.'s {r]t{k), k > 1} is well defined (which is not obvious in our setup) 
and constitutes an exchangeable sequence of {0, l}-valued random variables. We can 
then apply de Finetti's theorem, and prove that Xf^ — > Xt a.s for any fixed t > 0, 
and in probability locally uniformly in t, where Xt is a [0, 1]— valued Markov process, 
solution of the stochastic differential equation 



Xt = x-a f Xs{l-Xs)ds+ [ p{lu<x _ -X,-)M{ds,du,dp), (1.3) 

Jo ^[0,t]x]0,l[2 

where M{ds,du,dp) = M{ds,du,dp) — p^'^dsduA[dp), and M is a Poisson point mea- 
sure on M_|_x]0, l[x]0, 1] with intensity dsdup~^ K{dp) . Xt represents the proportion of 
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t 



Figure 1: The graphical representation of the A-lookdown model with selection of size N 
Solid lines represent type B individuals, while dotted lines represent type b individuals 



type b individuals at time t in the infinite size population. Note that uniqueness of a 



solution to (1.3) is proved in [8]. 

The paper is organized as follows. We both construct our process, and establish 
the crucial exchangeability property satisfied by the A-lookdown model with selection 
in section 2. In section 3 we establish the convergence of to the solution of (O). 
In section 4 we show that one of the two types fixates in finite time if and only if the 
A-coalescent comes down from infinity. Moreover, in the case of no fixation, we show 
that Xt — )■ Xoo G {0, 1} as i — )• oo, and give both a condition ensuring that X^o = a.s 
and a condition ensuring that F{Xoo = 1) > 0. Finally, we extend our results to the 
case A({0}) > in the last section 5. 

In this paper, we use N to denote the set of positive integers {1, 2, . . .}, and [n] to 
denote the set {1, . . . , n}. We suppose that the measure A fulfills the condition 

0< A((0,1)) <oo, A({1}) = 0, (1.4) 

and in all the paper except section 5, we assume that A({0}) = . 

2 The lookdown process, exchangeability 
2.1 Some results for general A 

Throughout the paper, the notation 

■■= / p'K{dp) 



is used for the rth moment of the finite measure A on [0, 1] for arbitrary real r. Note 
that fir is a decreasing function of r with oo > ^uq > > for r > 0, while fir 



may be either finite or infinite for r < 0. For r = 0, 1, ••• observe from (1.1) that 
fir = Xr+2,r+2 is the rate at which n„ jumps to its absorbing state {[n]} from any 
state with r + 2 blocks. Let X denote a random variable with distribution fi^^A, 
defined on some background probability space ($7, P) with expectation operator E, 
so E(X'' 



fir/fio- From (1.1), the transition rates of the A-coalescent are 



Xk,e = fioBiX'-^il - Xy-') for all 2 < £ < fc. 

For any partition with a finite number n > 2 of blocks, the total rate of transitions 
of all kinds in a A-coalescent equals 



n 



Xnj 



^ 1 — (1 — p)"' — np{l — p) 



n-l 



-A{dp) 



fioB 



l-(l-X)"-nX(l-X) 



71-1 



X2 



By monotone convergence. 



^[0,1] 



2.2 Construction of our process 

In this section, we will construct the process {rit{i),i > l,t > 0} corresponding to a 
given initial condition {r]Q{i),i > 1) defined in the Introduction. 



Recall the Poisson point process m defined in (1.2). For each n > 1 and t > , let 



I{n,t) = {k>l:tk£ [0,t] and #{It,,p, n [n]} > 2}. 

We have 

Lemma 2.1. For each n > 1 and t > 0, 

#I{n, t) < CO a.s 

Proof : Each atom {t,p) of m affects the n first individuals with probability 



Consequently 



The result follows. 



E(#/(n,t))< Qt^'A(dp) <oo. 



□ 
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2.2.1 A-lookdown model without selection 



For each > 1, one can define the vector = (^f (1), . . . , ^/^(A)), i > with values 
in {0,l}^by 

1. ^^(i) := r/o(i) for all i > 1. 

2. At any birth event {t,p) G m and such that {It,p Ci [A]} > 2, for each i G [A], 
Ct^{i) evolves as follows 

S,^ (min It^p), if i £ It,p \ {min It,p} 

C^- {i - {#{It,p n[l,...,i]}- 1)), otherwise. 

Using the above lemma, we see that the process has finitely many jump on [0, t] 

for all t > 0, hence its evolution is well defined. From this definition, one can easily 
deduce that, the evolution of the type at level i depends only upon the types at levels 
up to i. Consequently, if 1 < A < M, the restriction of to the first A levels yields 
in other words : 

{^f (1), . . . (A), i > 0} ^ {^f (1), . . .,^^{N),t > 0}. 

Hence, the process ry = is easily defined by a projective limit argument as a {0, 1}°°- 
valued process. 

2.2.2 A-lookdown model with selection 

This section is devoted to the construction of the infinite population lookdown model 
with selection. 

For each M > 1, we consider the process {'qf'^ {i), ^ > 1, i > 0) obtained by applying 
all the arrows between 1 < z < j < oo, and only the crosses on levels 1 to M. Using 
the fact that we have a finite number of crosses on any finite time interval, it is not 
hard to see that the process [r]^ ,t > 0) is well defined by applying the model without 
selection between two consecutive crosses, and applying the recipe described in the 
Introduction at a death time. More generally, our model is well defined if we suppress 
all the crosses above a curve which is bounded on any time interval [0, T]. Note also 
that, if we remove or modify the arrows and or the crosses above the evolution curve 
of a type B individual, this does not affect her evolution as well as that of those sitting 
below her. 

At any time t > 0, let Kt denote the lowest level occupied by a 5 individual. Of 

course, if Kq = l,Kt = 1, for all t > 0. If for any T, supo<t<T At < oo a.s, then the 
process {Kt,t > 0} is well defined by taking into account only those crosses below the 
curve Kt, and evolves as follows. When in state n> 1, Kt jumps to 

1. n + k at rate ("^^7^) A„+fc,fc+i, k > I; 

2. n — 1 at rate a{n — 1), a > 0. 
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In other words, the infinitesimal generator of the Markov process {Kt,t > 0} is given 
by: 

^9{'n) = '^i , jXn+k-i,k+i[9{n + k)-g{n)]+a{n-l)[g{n-l)-g{n)]. (2.1) 

fc=i ^ ^ 

Now, we are going to show that the process {r]t{i), i > 1, t > 0} is well defined. For 
this, we study two cases. 
Case 1: Kt — oo,t — oo. 
For each > 1, t > 0, we define 

= the level of the N— th individual of type B at time t. 

and 

= M{t > : i^f = oo}. 
We have > > • • • > 0. For each iV > 1, we define 

nN = {{s,ky,k<Kf}. 

Consider first the event 

^ = {ri^ = oo,Viv>i][l 

Recall the Poisson point measure m defined in ( |1.2[ ). Now, for each > 1, we define 
the process(r/^(i), i > l,t > 0), with values in {0, 1}°^, by 

1. t]q {i) '■= ??o(i) for all i > 1. 

2. At any birth event {t,p) G m,rj^^ evolves as follows 



rjf, (min h^p), if i £ It,p \ {min It^p} 

- i#{It,p n[l,...,i]}- 1)), otherwise, 



3. Suppose there is a cross on level j at time s. If ^ TIn or (s, j) £ Hn and 
r]g-{j) = 0, nothing happens. If G TIn and r]g~{j) = 1, then 



r]^-{i}, if i < j; 
r?f_(i + l), ifi>j. 



In other words, the process {ri^{i),i > l,t > 0} is obtained by applying all the arrows 



between 1 < i < j < oo, and only the crosses on levels 1 to Kj. . On the event A, we 
have a finite number of crosses on any finite time interval, and {r]^{i),i > l,t > 0) is 
constructed as explained above. Now, let 



^We shall see below that P(A) = 1 if the A-coalescent does not come down from infinity and P(A) = 
otherwise. 
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By a projective limit argument, we can easily deduce that the process {i]t{i),i > l,t> 
0} is well defined on the set T-L. Our model is defined on the event A. 

Now we consider the event A'^. We first work on the event {T^ < 00} . This means 
that the allele b fixates in finite time. It implies that for each N > 2, is finite as 
well. Consider first the process {r)l{i),i > l,t > 0} defined on Hi, i.e we take into 
account all the arrows between 1 < i < j < K}, and only the crosses on levels 1 to 
K}. This process is well defined on the time interval [0, T^^). However, on the interval 
[T;^, 00), r7j^(i) = l,Vi > 1, hence the process is well defined in "Hi. We next consider 
the process {r]f{i), « > 1, i > 0} defined on T-L2- This process is well defined on the time 
interval [0, T'^). But on the interval \T^,co), there is at most one B, whose position is 
completely specified from the previous step. Iterating that procedure, and using again 
a projective limit argument, we define the full A-lookdown model with selection. 

If = +00, but < +00 for some N, the construction is easily adapted to that 
case. In fact some arguments in section 4 below show that this cannot happen with 
positive probability. 

Case 2 : Kt 00, t ^ 00. 
Let 

Ti = inf{t >0: Kt = l}. 

It is not hard to deduce from the strong Markov property of the process Kf that 
{Ti < 00} a.s on the set {K^ 00, t — )• 00}. In that case the idea is to show that 
there exists an increasing mapping ijj : N ^ N such that a.s. for N large enough, 
any individual sitting on level ^{N) at any time never visits a level below N, with the 
convention that if that individual dies, we replace him by his neighbor below. Once 
this is true, the evolution of the individuals sitting on levels 1, 2, . . . , A'' is not affected 
by deleting the crosses above level ip{N). Hence it is well defined. If this holds for all 
N large enough, the whole model is well defined. 
Let 

M = sup Kf. 

0<t<Ti 

For each N > M, let ip{N) = Ne""^^ {Ne"^^"^ + 1) + Kq, where Sn is defined below. 
We want to show that an individual sitting on level (p{N) + at any time t > never 
visits a level below N. In order to prove this, we couple our model with the following 
one. 

On the interval [0,Ti], we erase all the arrows pointing to levels above Kt, and 
pretend that all individuals above level Kt, < t < Ti, are of type b, i.e coded by 1, 
and we apply all the crosses above level Kt. This model is clearly well defined since 
until Ti there is only one 0, all other sites being occupied by I's. We next extend this 
model for t > Ti as follows : 

For each t >Ti, let Kt denote the lowest level occupied by a h individual. At time 
Ti, ?7ri(l) = 0, T]Ti{i) = 1, for all i > 2. At any time t > Ti, we shall have T]t{i) = for 
i < Kt, and r/t(i) = 1 for i > Kf. Again all crosses are kept, and we keep only those 
arrows whose tip hits a level j < Kf 

This model is well defined. For each > 1, wc define by Sn the first time where 
all the N first individuals of this model are of type B. We have 
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Lemma 2.2. IfTi < oo, then for each N >1, 

Sn < oo a.s 

Proof : The result follows from Ti < oo and the fact that the process of arrows from 
1 to 2 is a Poisson process with rate A((0, 1)). 

□ 

Now, let t > 0} denote the process which describes the position at time t 

of the individual sitting on level ip{N) at time in the present model. 

The individual who sits on level (p{N) at time will remain below the level (p(N)+N 
on the time interval [0, S'Ar]. If she does not visit any level below N before time Sn, 
she will never visit any level below N at any time, and moreover any individual who 
visits level <p{N) + before time will remain above the individual who was sitting 
at level ^p{N) at time until Sn, hence will never visit any level below A^. 

Since the "true" model has more arrows and less "active crosses" than the present 
model, if we show that in the present model a.s. there exists N such that the individual 
who starts from level 93 (A^) at time never visits a level below N, we will have that in 
the true model a.s. for N large enough the evolution within the box {t, i) G [0, 00) x 
{1,2,..., N} is not altered by removing all the crosses above '^{N) + N. A projective 
limiting argument allows us then to conclude that the full model is well defined. 

The result will follow from the Borel-Cantelli and the next lemma. 

Lemma 2.3. IfTi < 00, then for each N >M, 

Pjv(30 < i < 5iv such that ^f^^ < N) < j^, 
where PnI]= P(- | Sn) 

Proof : It is clear from the definition of ^f^^^ that there exists a death process 
{Dt,t > 0), which is independent of {Kt,t > 0) conditionally upon Dq = ip{N) — Kq, 
and such that 

^fW=^, + A, \/t>o, 



where 



< t < Ti; 




On the other hand, we have 



Mm 

All we need to prove is that 



^n<^^. ^* > ^> ^ ^n<^^. A > A^} D {Ds^ > N}. 
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The process {Dt, t > 0) is a jump Markov death process which takes values in the space 
{0, 1, . . . , ^{N) — Kq}. When in state n, Dt jumps to n — 1 at rate an. In other words 
the infinitesimal generator of {Dt,t > 0} is given by 

Qf{n) = an[f{n-l)-f{n)]. 
Let / : N — )• M. the process {Mf)t>o given by 

m/ = /(A) - f{Do) -a f DsifiDs - 1) - f{Ds)]ds (2.2) 

Jo 



is a martingale. Applying (2.2) with the particular choice /(n) = n, there exists a 
martingale (M/)t>o such that Mq = and 

Dt = Do-af Dsds + Ml, t > 0. (2.3) 
Jo 

We note that {Ml,t > 0} is a martingale under PAr[.]. This is due to the fact that the 
Poisson process of crosses above Kt is independent of Kt. We first deduce from (2.3) 
that Bn{Ds) = Doe-"'''. 

Using the fact that Dt is a pure death process, we obtain the identity 

[M% = Do- Dt, 



which, together with (2.3), implies 

< >t= a [ 
Jo 



Dgds. 



From (2.3), it is easy to deduce that 



Dt = e-'^\^{N) -Kq)+ [ e-"(*-^)(iMi, 

^0 

which implies that 



PTvfl / e'^'dMll > Ar2e2°-5iv 







1 f^^ ^ 

2 



The result is proved 



□ 



From now on, we equip the probability space (il,J-", P) with the filtration defined 
by J't = o'{r/s(i), i > 1, < s < t}. Any stopping time will be defined with respect to 
that filtration. 
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2.3 Exchangeability 

In this subsection, we will show that the A-lookdown model with selection preserves 
the exchangeability property, by an argument similar to that which we developed in 

m- 

Let Sn denote the group of permutations of the set {1,2, ... ,n}. For all it £ Sn 
and a["l = {ai)i<i<n G {0, 1}", we define the vectors 

7r~^(a["]) = (a^-i(i), . . . ,a^-i(„)) = (af )i<i<„ 

^(el"b = (6Kl)),...,6(vr(n)) 

We should point out that vr(^|"') is a permutation of (■^t(l), • • • , ^ti^)) and it is clear 
from the definitions that 

{7r(el"b=a'"l} = {ef'=vr-\aW)}, for any n e Sn- (2.4) 

The main result of this subsection is 

Theorem 2.4. If (?/o(0)i>i ^''^^ exchangeable random variables, then for all t > 0, 
{rjt{i))i>i are exchangeable. 

We first establish two lemmas 

Lemma 2.5. For any finite stopping time r, any ^-valued Tr-measurable random 
variable n* , if the random vector r^r ' = (?/t(1)i • • • i Vrin*)) is exchangeable, and T is 
the first time after t of an arrow pointing to a level < n* or a death at a level < n* , then 
conditionally upon the fact that T = ti^, for some io>l and Itig,pig H [n*] = k, where 

k > 2, the random vector t?!"^ ^^^''^ = [rjn^ (1), . . . , rjt^^ (n* - 2 + /c), r/t-^ (n* - 1 + k)^ 
is exchangeable. 

Proof : For the sake of simplifying the notations, we condition upon n* = n, ti^ = t, 
PUg = P a-iid Iti^,pi^ n [n*] = k. We start with some notation. 

^jovjfc-i ._ i^j^g ^ levels selected by the point {t,p) between levels 1 and n are 
jo,Ji, ■ • • ,jfe-l}- 

We define 

Vt,n[] =P(-|iio =t,n* = n,It,pr\[n] = k). 
Thanks to we deduce that, for vr G Sn-i+k, a^'^-^+^'^ G {0, 



l<io<ii<--<jfe~i<n 

On the event A^''''"'^'^'^ , we have : 



Vt-{i), if 1 < « < Ji 

rit-Uo), if i G {ii,j2, • • -Jk^i} 

^m-{i - (#{{J1,J2, • • • , jfc-i} n [i]})), if ji <i<n-l + k, i^ {j2, . . . , jfc-i}. 
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This implies that 



For 1 < ii < J2 < • • • < jk^i < n, define the mapping Pj^,j2,-,jk-i '■ {0' 1}"+''"^ 
{0, 1}" by : 

Pjij2,.-Jfc-i(^l' • • • ' ^n-l+fc) = (^jii • • • 1 ^jfe-i)) 

where 



^im = • • • , - 1 ) , 2 < m < fc - 1 

= (ftjfe-i+l) ^J*-i+2, • • • , bn-l+k)- 
In other words, Pji,j2,--,jk-ii'^) vector z from which the coordinates with indices 



ji, . . . ,jk-i have been suppressed. The right hand side of (2.5) is equal to 



i<io<ii<--<ifc-i<™ 

It is easy to see that the events {rj-^} = yOjij2v -ifc-i(^^^('^'" "'^^'^'))) ^'^^ _^o,---,ik-i 
independent. Thus 



^ 30 n 3k-l' 



i<io<ii<--<ifc-i<n 
n\ 

, ] / ^ Ija'^ =a'^ ■■■=a'^ | 

A; / '^-^ n n ik-i' 

i<io<ii<--<ifc-i<™ 

On the other hand, we have 

#{1 < io < • • • < jk^i < n : = • • • = aj^_-^} = #{1 < jo < • • • < jk-i < n : a]^ ■ ■ ■ = aj^_ 

Let io < ii < ■ ■ ■ < ik-i be the increasing reordering of the set {vr(jo), vr(ji), • • • , 7r(jfc_i)}. 
If = = • • • = aj,_,, then we have aj^ = aj^--- = a^^_^ = aj^ = aj^ = ■ ■ ■ = 

a 



^, and consequently Pji,j2v jt-i (^^" i+'^l) and Pei,e2,-/k-ii'^ ^(a'" i+'^l)) contain 
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\n\ 

the same number of O's and I's. Since r/^J is exchangeable, 

^ ^ 7e{0,l} l<^Q<^i<-<4-i<n 
^ ^ 7G{0,1} l<JO<Jl<-<Jfc-l<n 

The result follows. □ 

Lemma 2.6. For any finite stopping time t, any l>i -valued -measurable random 

variable n* , if the random vector rj^ ' = (r/T-(l), . . . , r/^(n*)) is exchangeable, and T 
is the first time after r of an arrow pointing to a level < n* or a death at a level 
< n* , then conditionally upon the fact that T is the time of a death, the random vector 
Vt ~ {vt{^), ■ ■ ■ 1 ririn* — 1)) is exchangeable. 

Proof : To ease the notation we will condition upon n* = n and T = t. Let vr G Sn-i 
be arbitrary. We consider the events : 

Bl := {the level of the dying individual at time i is i}. 

Let Pt,„[.] = P(.|T = t,n* = n). We have 

Pt,n(7r(r?r'l) = at-^l) = E h,n [rlr'^ = ir-\a^--% Bl) 

l<i<n 

= E P(^t(l) = a^...,r/t(n-l) = <_i,S|). 

l<i<n 

Define 

c^" = (a^,...,a^_i,l,af,...,<_i), = (ai, . . . , Oj-i, 1, a^, . . . , a„_i). 

The last term in the previous relation is equal to 

E ^t,n (^J = Bl) = E p (^-' = cr) hn {bi i ^fJ = <-) 

l<i<n l<i<n 

Thanks to the exchangeability of (r/^- (1), . . . , (n)), we have 

J- + Z^i=i «i i<i<„ 
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since l^J=i ci] = Yl^jLi^ ^^'^ ^ permutation of cj. The result follows. □ 

We can now proceed with the 



Proof of Theorem 2.4 For each > 1, let {V/^ , ^ > 0} denote the N-valued 
process which describes the position at time t of the individual sitting on level at 
time 0, with the convention that, if that individual dies, we replace him by his neighbor 
below. The construction of our process {rjt{i),i > l,t > 0} in section 2.2 shows that 
t>0 Vt 

that for each t > 0, iV > 1, (ryt(l), . . . , rit{Vt^)) 



It follows from Lemma 



2.5 



and 



2.6 



is an exchangeable random vector. 

Consequently, for any t > 0, n > 1, vr G 5„, a'"! G {0, 1}, 

|P(r?!"] = aN) - P(,,M = 7r-i(aW))| < P(y,^ < n) 

which goes to zero, as — ?• oo. The result follows. 



□ 



For each A'^ > 1 and t > 0, denote by the proportion of type b individuals at 
time t among the first A'^ individuals, i.e. 

1 ^ 

^^ = nY.^^^^- (2-6) 

i=l 

We are interested in the limit of {X^^)t>o as A^ tends to infinity. For this, let us recall 
the following useful result due to de Finetti (see e. g. [2]). In this statement, Q denotes 
the tail cr— field of the sequence {Xn,n > 1}. 

Theorem 2.7. An exchangeable (countably infinite) sequence {Xn,n > 1} of random 
variables is a mixture of i.i.d. r.v.'s, in the sense that conditionally upon Q, the Xn's 
are i.i.d. 

As a consequence, we have the following asymptotic property for fixed t of the 
sequence (X^)n>i defined by Q 



Corollary 2.8. For each t >0, 



Xt = lim exist a.s. (2.7) 

3 Tightness and Convergence to the A-W-F SDE 
with selection 

3.1 Tightness of {X^, iV > 1, t > 0} 

In this part, we will prove tightness of {X^)j~]>i in D([0,oo[), where for each A^ > 1 
and t > 0, is defined by K^. We start with some notations. 
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For any N, n, r,p such that N > l,Nr eN,r g]0, G [0, 1], we define 
Y{-, N,p) to be the binomial distribution function witli parameter N andp; H{-, N, n, r) 
the hypergeometric distribution function with parameter (A^— l,n— 1,^^^); H(;,N, n, r) 
the hypergeometric distribution function with parameter ( A'' — 1 , n — 1 , -^^j ) • For every 
v,w e [0, 1], let 

mi{s:Y{s.N,p) > v}, 
inf{s; H{s, N, n, r) > w}, 
inf {s; H{s, N, n, r) > w}. 

It follows that if V,W are W([0, 1]) r.v.'s, then the law of F^{V) is binomial with 
parameter N,p. G N.n.riW) {lesp G]\f,n,r(W)) is hypergeometric with parameters A'" — 
(resp Af-l,n-l,^'). Note that = Y~H; N,p),GN,nA-) = 

H~^{-, N,n,r) and GN,n,r{') = H~^{-, N,n,r). We recall that if X is hypergeometric 
with parameters {N, n,p) such that Np G N and p G [0, 1], then 

N — n 

E(X) = np and Var{X) = — — -np{l —p). 
Now, for every r, u,p, v,w & [0, 1], let 

Prom the identity r(n — 1 — 'E[GN,n,r{w)]) = (1 — r)'E[GN,n,r{w)], we deduce the 
Lemma 3.1. For each N > l,r,p,v G [0, 1] and t>0, 

/ tp^ {r,u,p,v,w)dudw = 0. 



GN,nAw) 



i) {r,u,p,v,w) = —1^^ 



Using the definition of the model, it is easy to see that 

X^ = X^+ [ il;^{X^-,u,p,v,w)Mo{ds,du,dp,dv,dw) 

i[0,t]x]0,l]4 



-'^ J[0,t]x[0,l]2 - - s- 

where Mq and Mf^ are two mutually independent Poisson point processes. Mq is 
a Poisson point process on M_|_ x [0, 1] x [0, 1] x [0, 1] x [0, 1] with intensity measure 
IJ,{ds,du,dp,dv,dw) = dsdup^^A{dp)dvdw, is a Poisson point process on 1R+ x 
[0, 1] X [0, 1] with intensity measure aNdsdudv. Now, let 

Mo = Mo - n, = Mf - aN\. (3.1) 
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Using Lemma 3.1, we have 



[0,t]x]0,l]'' 



Tp (X^- , f , w)MQ{ds, du, dp, dv, dw) 



N 



K<x^_ l^>x^+i^f (^«' du, dv)-a I Xf (1 - X;^+')ds 



[0,t]x[0,l]2 



(3.2) 



For each N > 1, t > 0, we define 



TWf = / ip^{X^^,u,p,v,w)Mo{ds,du,dp,dv,dw) 



[0,t]x]0,l]* 



1 



Mr = — l„<^^ l^>^iv+iMf {ds, du, dv) 

i[0,t]x]0,l]2 - - 



v!' = -a fxf{l-X^^')ds. 
Jo 



Aif and AAj are two orthogonal martingales. It is easy to see that 

X[^ = X,^ + V,^ + M^-Mr. 

\/N > i,X^ E [0, 1], which imphes that it is tight. Moreover, we have 
Proposition 3.2. The sequence {X^,N > 1) is tight in D{[0,oo]). 

We first establish the lemma : 
Lemma 3.3. For each N > I and t > 0, 



{Mnt = Ami)) / Xf(l-Xf)ds 
Jo 

ft 



{M% = ^I^X^il-X^^')ds 

Proof : Using the fact that and are pure-jump martingales, we deduce that 

{xp'^ (X^ , u,p, V, w))"^ dsdup"'^ A{dp)dvdw . 



{M% 



[0,t]x]0,l]'' 



Let 



A^{X^,p)= [ {ij^{X^,u,p,v,w))^dudvdw 
J]o,ip 

1 



X^ (Fp^'iv) - 1 - G^,i.iV(„),;,iV (w) 



+ {l-X^){G^^pN^,^^xi^iw)y 



dvdw. 
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Tedious but standard calculations yield 



[0,1] 



riF^^{v) - 1 - G^^pN(^^)^riw)j dw + {l- r) f ) 



^r(l-r)F,^(r;)(F,^(r;) 



for every v,r £ [0, 1]^. Consequently 
•A (A, ,p)- 



[0,1] 



We deduce that 



{M% = [ A''{X^,p)dsp-^k{dp) 

J[0,t]x[0,l] 

= A((0,1)) / Xf(l-Xf)ds. 
J [OA 



Similarly, we have 



a 
N 
a 
N 



[0,t]x[0,l]x[0,l] 



/ (1 - 

J [OA 



The lemma has been established. 

We can now proceed with the 
Proof of Proposition 13.21 Let 



^{x,y) = -ax{l - y) 



a 



^N{x,y)= ^A((0,l))(l-x) + -(l-y) 



It follows from lemma 3.3 that for each t > 0, 

rt 



Xf = Xo+/ ^{X^,xr')ds + Mf -Mi' 
Jo 



and 



Jo 

IVloreover, VT > 

sup sup(|(^(Xf,Xf+i)|+*^(Xf,Xf+i))<C7 a.s. 

0<t<T N>1 ^ ' 
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Aldous' tightness criterion (see Aldous |1J) is an easy consequence of (3.3). □ 



Now, from Proposition 3.2 and (2.7), it is not hard to show there exists a process 



X G D{[Q,oo)), such that for ah t > 0, 

Xl^^Xt a.s, (3.4) 

and 

X^ ^ X weakly in D{[0, oo)). 
We deduce for the above facts 

Theorem 3.4. For all T > 0, 

sup \X[^ — Xt| — )• in probability, as N ^ oo. 

0<t<T 

Proof : To each 5 > 0, we associate n > 1 and = to < *i < • • • < = 7", such that 
supi<j<„(ti — ti-i) < 6. We have, with the notation y A z = inf(y, z), 

sup |Xf -Xt\< sup sup |Xf - J A |Xf -Xi^\+ sup - | 

0<t<T i ti-i<t<ti i 

+ sup sup \Xt - Xt^\ A\Xt - Xt._^ I 



< wUx'^ , 6) + sup - XiJ + «;^^(X, 6), 



where 



w''^{x,6)= sup |x(t) — x(ti)| A |2;(t) — x(t2)|- 

0<ti<t<t2<T,t2-ti<S 



Since from (3.4), 



sup \X[^ — Xt^ \ — )• a. s., as — )• oo, 

l<j<n 

limsupPf sup iXf-Xil > e) < hmsupPfw^(^^> > e/2)+p(wT{X,6) > e/2 

From the proof of Theorem 13.5 in fT], we know that tightness of imphes that 

limlimsupP(u;^(X^,(5) > e/2) = 0. 

Since X is cadlag a.s., for each e > 0, 

P{wt{X, S) > e/2) ^0, as (5 ^ 0. 
The result follows. □ 



3.2 Convergence to the A- Wright-Fisher SDE with selec- 
tion 



Our goal is to get a representation of the process Xt defined in (3.4) as the unique 
weak solution of the stochastic differential equation (1.3). 

Let (ri, T, P) be a fixed probability space with a filtration {Tt)t>o satisfying the 
usual conditions. Recall the Poisson point measure M = ^ti,Ui,pi defined in the 

Introduction, and for every u g]0, 1[ and r £ [0, 1], we introduce the elementary function 

^'('U,, r) = lu<r - r. 



We rewrite equation (1.3) as 



Xt = x-a / Xs{l-Xs)ds + 



[0,t]x]0,l[2 



{u, Xs^)M {ds, du, dp) , t > 0,0 < X < 1, 

(3.5) 

which we call the A- Wright-Fisher SDE with selection. Without loss of generality, 
we shall assume that a > 0, which means that Xt represents the proportion of non- 
advantageous alleles. 

The proof of the following identity is standard and left to the reader. 

Lemma 3.5. For each r £ [0, 1], 



[0,1]* 



[ip {r,u,p,v,w) — p"^{u,r)^ dup A{dp)dvdw 



2r(l - r) 



N 



IN - 1 



(1 - up)^-^duK{dp) 



A([0,1]) 
A^- 1 



We are now in position to prove our main result. 



Theorem 3.6. Suppose that Xq — )• x a.s., as N ^ oo. Then the [0, 1] — valued process 
{Xt,t > 0} defined by (3.4) is the (unique in law) solution of the A-W-F SDE with 
selection (3.5). 



Proof : Uniqueness of the solution of (3.5) follows from Theorem 4.1 in [8J. We now 



prove that {Xt)t>o defined by (3.4) is a solution of the A- Wright-Fisher (3.5). 

7-7 



Since Xj^ — )• Xt a.s., from (|3.1|), it suffices to show that for each i > 0, as N ^ oo 



Vt^ -Mt^ ^-a / X,{l-X, 







s)ds+ j 
Jk 



[0,t]x[0,l]2 



p^'(X^- , u)M{ds, du, dp) 



in probability. Now using (3.4) and Lemma 3.3, it is easy to see that, for each t>0 

Vt^ —J' -a Xs{l - Xs)ds a.s, and Aff — )• in probability, 
Jo 

as A" — )• oo. The main work is to show that 



M 



N 



[0,i]x]0,l[x]0,l[ 



p'${Xt- , u)M{ds, du, dp) in probability, as N ^ oo. (3.6) 
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For each X > 1 and t > 0, let 

h^it)= [ {ip^{Xf.,u,p,v,w) -p^{X,-,u)) Mo{ds,du,dp,dv,dw), 



where Mq is defined by (3.1). h {t) is a martingale, and 

{h^)t = [ {il)^{X^ ,u,p,v,w)-p^{Xs,u)f dsdup-^K{dp)dvdw. 

J[0,t]x[0,l]4 



We have 

{h^)t<2 [ {p^iu,X^)-p^{u,Xs))^dsdup-^A{dp) 

J[0,t]x[0,l]2 



+ 2t sup / {'ip^{X^,u,p,v,w) -p^{u,X^))'^ dup-^A{dp)dvdw 
o<s<t y[o,i]4 

<2 / {p'^{u,X^) -p^{u,Xs))^ dsdup~^A{dp) 

J[0,t]x[0,l]2 



+ 2t sup / V, ti)) — r)) dup '^A{dp)dvdw. 

0<r<l J[0,1]'' 

Using the fact that X^ — )• Xs a.s., it is not hard to show by the dominated convergence 
theorem that as X — )• oo, 

/ p^{^{u,X^)--^{u,Xs)f dsduA{dp)^{) a.s. (3.7) 

J[0,t]x[0,l]2 

Now from lemma 3.5 , it is easy to show that as X — t- oo, 

sup / [il^^ {r,u,p,v,w) — p"^{u,r)Y dup~^A{dp)dvdw ^ (3.8) 

0<r<l J[0,l]4 



Combining (3.7) and (3.8), we deduce that 

yt > 0, {h^)t a.s, as N -^oo. 
On the other hand, we have 

{h^)t < CtA([0,l]), VX > 2. 
We deduce from the dominated convergence theorem that 

lim E[h'^{tf] = Vt > 

i.e 

Mf = f ^^{X^.,u,p,v,w)Moids,du,dp,dv,dw) ^ [ p^iX,-,u)M{ds,du,dp) 

J[0,t]x[0,l]4 J[0,t]x[0,l]-^ 
as X —7- oo, in particularly 

—7- / p'if{Xg- , u)M{ds, du, dp) in probability , as N ^ oo . 

i[0,t]x[0,l]2 



(3.6) is established. 

□ 
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Remark 3.7. Uniqueness in law could also by proved as in l^], (where the case q = 
is treated) by a duality argument, which we now sketch . 

Recall the notation ^{u,y) = lu<y — U- For every y G [0,1] and every function 
5 : [0, 1] — )• M of class C^, we set 



[0,1] X [0,1] 



g{y+p^{u,y))-g{y)-p^{u,y)g'{ii) p '^A{dp)du-ag'{y){l-y)y. 



A solution {Yt)t>o of (3.5) is a Markov process with generator C Hence for every 
5 : [0, 1] — )• M of class C'^, the process 



g{Yt)- f dsCg{Ys), t>0 



is a martingale. 

It is plain that for g{z) = z"" 



n—k+l n 



(3.9) 



to 



Let {Rt,t > 0} be a N -valued jump Markov process which, when in state k, jumps 

1. k-l+l at rate {f)\k,e, 2<i<k; 

2. k + 1 at rate ak, a > 0. 

In other words, the infinitesimal generator of {Rt,t > 0} is given by: 

J^*f{k) = Q XkAfik -£+!)- f{k)] + ak[f{k + 1) - f{k)]. 



For every z £ [0, 1] and every r € N, we set 

P{z,r) = z\ 
Viewing P{z,r) as a function ofr, we have 



(3.10) 



k=2 ^ ^ 

On the other hand, viewing P{z, r) as a function of z we can easily evaluate CP{z, r) 
from formula (3.9), and we deduce that 

CP{z,r) = C*P{z,r). (3.11) 



Now suppose that {Yt)t>o is a solution of (3.5), and let Rq = n. By standard argument 
(see Section 4-4 deduce from (3.11) that 

np{yt,R^)] = np{Yo.Rt)], 
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I.e 



E[y-|yo 







n 



Since this is true for each n > 1 and Yt take values in the compact set [0, 1], this is 
enough to identify the conditional law of Yt, given that Yq = x, for all < x < 1. 
Since (lt)t>o 'is a homogeneous Markov process, this implies that the law of (lt)j>o is 
uniquely determined. 



If we prove a priori that {Xt)t>o defined by (3.4) is a Markov process, we can use 



the fohowing Remark to prove that Xt is a weak sense solution of the A- Wright-Fisher 



SDE (3.5). 



Remark 3.8. Suppose we know that Xt defined by (3.4) is a Markov process. Let us 



look backwards from time t to time 0. For each < s < t, we denote by Z"'* the highest 
level occupied by the ancestors at time s of the n first individuals at time t. We know 
that conditionally upon Xt, the {r]t{i),i > 1} are i.i.d Bernoulli with parameter Xt. 
Consequently, for any n > 1, 

Xr = P(r?i(l) = ...=r/i(n) = l|Xi), 

this implies that 

E,[Xf] = B,\P{rjt{l) = ■■■ = r,t{n) = 1 | Xt)] 
= P,(r?t(l) = --- = r?t(n) = l) 
= 'Px{the 1 . . . Zq'* individuals at time are all b) 

rrn.t 



It is plain that the conditional law of Zq' , given that {rit{l) = ■ ■ ■ = r]t{n) 
the conditional law of Rt, given that Rq = n. Consequently, for each n > 1 



1) equal 



E[Xr \Xo = x] 



Ya 



where {Yt)t>o is a solution of (3.5). But for all t > 0,r € [0, 1], the conditional law of 
Xt, given that Xq = x is determined by its moments, since Xt is a bounded r. v. So 
Xt and Yt have the same transition densities, that is {Xt,t > 0} is the unique weak 



solution of (3.5) 



4 Fixation and non-fixation in the A-W-F SDE 
4.1 The GDI property of the A-coalescent 

In this subsection, we recah a remarkable property of the A-coalescent (Ilt)t>o defined in 

In] ~ 

the introduction. For each n > 1, let #11^ denote the number of blocks in the partition 
jjM (jjW jg ^j^g restriction of Ut to [n] ). Then let r„ = inf{i > : #n["' = 1}. As 
stated in (31) of [T5], we have 

= Ti < T2 < T3 < . . . t Too < oo. 
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We say the A-coalescent comes down from infinity (A G GDI) if P{^Ilt < oo) = 1 
for all t > 0, and we say it stays infinite if P(#nt = oo) = 1 for all t > 0. In terms 
of the population model, this means that for any i > 0, we can find a finite number 
of individuals in the initial population which generate the entire population at time 
t. The coalescent comes down from infinity if and only if Too < oo a.s. We will show 
that this is equivalent to fixation. Kingman showed that the (5o-coalescent comes down 
from infinity. 

A necessary and sufficient condition for a A-coalescent to come down from infinity 
was given by Schweinsberg [T7]: define 



4>{n) 



k=2 ^ ' 



and 



u{dp) = p ^A(dp). 
It is not hard to deduce from the binomial formula 



(j){n) 



[ [np - 1 + (1 -pT]iy{dp). 
Jo 



Schweinsberg's result |17J says that the A-coalescent comes down from infinity if and 
only if 



n=2 



< oo. 



(4.1) 



The condition of convergence of this series is also necessary and sufficient for fixation 
in finite time. Using the fact that the function fn{p) = (1 — p)^ — 1 is decreasing for 
any fixed n, we have 

/■I 

{np — l)v{dp) < (t){n) < n pv{dp), Vn > 1. 

Jo 

The last assertion together with (4.1), implies that if Jq ph'{dp) < oo then the A- 
coalescent stays infinite. The result has been proved by Pitman (see lemma 25 in 

m)- 

Theorem 3.6 shows that Xt is a bounded supermartingale. Indeed, if {Xt)t>o is a 

p^{u,Xg-)M{ds,du,dp) I 

[0,t]x]0,l[x]0,l[ 



solution of 3.5 , then for all < t < s. 



E{Xt\Ts)<x-a Xr{l- Xr)dr + B 
Jo 

= Xs. 

Consequently the following limit exists a.s 

X, 

Indeed, and 1 are the only possible limit values. 



oo - lim Xt e {0, 1}. 



(4.2) 
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4.2 Fixation and non-fixation in the A-W-F SDE 

We assume that the initial proportion x of type B individuals satisfies < x < 1. In 
this section, we prove that if the condition (4.1) is satisfied, then we have fixation in 
our model in finite time. Before establishing the main result of this section, we collect 
some results which will be required for its proof 



Lemma 4.1. 

where v{dp) =p^^A{dp). 
Proof : 

4>{n) ■■ 



(bin) 

t / P^{dp) as n t oo, 



n 



[ [np-l + (l-p)"]i^(dp) 
Jo 

(1 - up)''-^du 



A{dp). 



On the last line, we have made use of the identity 

(l-p)"-l= / -np{l-up)''-^du. 
Jo 

For each p g]0, 1], let 

1 



rip) = -11 

P 



(1 - up^-'du 



We have. 



n-V(n) = / r{p)Hdp). 



The result follows from the monotone convergence theorem. 

We now deduce that 
Lemma 4.2. The function (p increases, and 

oo ^ 

— — < oo =^ — - 

fOlni ^ — ' (Bin 

n=2 



oo ^ 

0(n) 



(n) — an 



< oo. 



Proof : We have 



(n+ 1) - <f>{n) 



Jo 



[i-pTHdp) 



p{i-{i-pr)v{dp) 



□ 



> 0. 
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Which imphes the first claim. Now, we already know that if X]jS=2 ^ then 

pu{dp) = OO. Thus, the second assertion is a consequence of the last lemma and the 
following relation 

CXD ^ OO OO ^ 

J- v-^ Ci s—^ i 

^-^ (j){n) — an ^-^ (p{n){n~^(j){n) — a) ^-^ ^ij>) 



□ 



The lemma is proved. 

For each t > 0, we define again 

Kt = inf{i > 1 : r,t{i) = 0}. 

and 

Ti = inf{t >0: Kt = l}. 

We have the following 

Theorem 4.3. If A & CDI, then one of the two types ( b or B) fixates in finite time, 
i.e. 

3 C < OO a.s : = Xoo G {0, 1} 

If Ki CDI, then 

Vt > 0, < < 1 a.s. 

Proof : The proof has been inspired by [S] (see Section 4 ). 

Step 1 : Suppose that A G CDI. We consider two cases. 
Case 1: Kq = 1. 

In this case, the allele B fixates in the population. Indeed, the individual at level 1 
never dies and, he cannot be pushed to an upper level. Let 

C = inf{t > : r}t{i) = 0,Vi > 1}. 

C is the time of fixation of allele B. We are going to show that C < oo a.s. 

Let > 1 denote a fixed integer, and we define Sn the first time when the first 
individuals are all of type B. As Sm < oo a.s, for all > 1, we wait until all the first 
individuals are of B type and we look backward in time, and let Z^' denote the number 
of ancestors at time t of the A^ first individuals at time (i.e of the individuals sitting 
on levels 1, ... ,A^ at time Sn)- The process {Zl^)t>o is a jump Markov process with 
state-space {1,2,..., A^}. When = n, Z^ is shifted to n — £ + 1, where 2 < £ < n, 
at rate (")A„^£. In other words, the infinitesimal generator of Z^ is given by 

^'^^(") = E (") ^«.^[/(^ - ^ + 1) - /HI- 

Let 

Ctv = inf{t > : = 1}. 
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Ctv is also an upper bound for the time taken in the forward time direction, for the 
progeny of a type B individual sitting on level 1, to invade all the levels {1,2,..., N}. 
The sequence {Cn)n>i increases and it is easy to see that limAr_>oo Cn = C, where ( is 
the time of fixation of the type B. We already know that since A G GDI, 

V — < oo. (4.3) 

Now, for each n > 1, we define 

oo ^ 
k=n+l ' 

We have for 2 < ^ < n 

/(„-f+i)-/(„)= ± J^. 

k=n-l+2 ' 



Recall lemma 4.2 Since Xj^ is decreasing, we have for 2 < £ < n, 

/(n-^ + l)-/(n)>(^-l)- ^ 



and therefore 



(£-l)A„,£ = l. 



Using the fact that the process 

/(Zf ) - /(iV) - r Q^/(Zf t > 

Jo 

is a martingale, we obtain 

= /(l)-/(A^) 
</(l) 

The last assertion, together with the monotone convergence theorem, implies that 

E(C) < /(I) < oo. 



Case 2 : Ko> 1. 

If Ti < oo then B fixates in finite time. Indeed, wait until Ti which is a stopping time 
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at which the Markov process {r]t{i),i > l}t>o starts afresh, and then use the argument 
from Case 1. 

We suppose now that Ti = oo, which impUes that Kt — >■ oo,t — >■ oo, as already 
noted in section 2.2.2. In other word, if Ti = oo, then the ahele B does not fixate in 
the population. Let 

no = inf{n > 1 : (j){n) — an > 1}. 

Such an no exists because since A G GDI, pi'{dp) = +oo, hence by Lemma 4.1, we 
have lim„_j.oo n~^(j){n) = +oo. Let N > n^ denote a fixed integer. We let 

5^ = mi{t > : rjtil) = ■■■ = miN) = 1} 

Since Kt — )• oo, t — )• oo, we have that Sn < oo a.s. Now, we are going to use a similar 
argument as in Case 1 . We wait until the N first individuals are all of b type and 
we look backward in time, and we denote by the highest level occupied by the 
ancestors at time t of the N first individuals at time Sn- The process {Zl^)t>o is a 
jump Markov process with state-space {l,2,...,oo}. When Zf = n, Zf jumps to 
n — £ + 1, where 2 < ^ < n, at rate (^)Xn,e, and jumps to n + 1 at rate an. The 
infinitesimal generator of is given by 

^'^^(^) = E - ^ + 1) - f{n)] + an{f{n + 1) - /(n)). 

Let 

C = inf{t > : Zf < no}. 
Now, we are going to show that the allele b fixates in finite time. For this, we need 
only prove that ( = lim7v-s>oo Cn" < °° 

Recall lemma |4?2| For each n > 1, we define 

oo ^ 

" ^ U(k) - ak) V l 

k=n+l ' ' 



By Lemma 4.2, for each n > 2, /(n) is finite. We have for 2 < £ <n 

1 



f{n-£ + l)-f{n)= 



Since k — )• l/{(p{k) — ak) V 1 is decreasing, we obtain 

f{n-£+l)- f{n) >{£- I) ^ 

[(pyn) — an) V 1 

and therefore 

{4>{n) -an)ylj-i \£J 



=2 

4>{n) an 



(n+ 1) - a(n + 1)) V 1 



((/)(n) - an) V 1 {(t>{n + 1) - a{n + 1)) \/ 1 
^ 4>{n) an 

" {4>{n) — an) V 1 — ctn) V 1 
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hence f{n) > 1, for each n > no- Using the fact that the process 

/(Zf ) - /(AT) - /* Q^/(Zf t > 

Jo 

is a Martingale, we obtain 



= /(Zf„o)-/(iV) 

</(l) 

Since the sequences C]^ t C as t , the last assertions, together with the monotone 
convergence theorem, implies that 

E(C) < /(I) < oo, 

and hence ( < oo a.s. 

Step 2 : Suppose A ^ GDI, that is the A-coalescent does not come down from 
infinity. We have 



n=2 '^(") 



+00. (4.4) 



We claim that {Kt,t > 0) does not reach oo in finite time. The contrary would imply 
that 3T < oo such that Kt = oo a.s., so the number of ancestors of the infinite 
population at time T of the A-lookdown model would be finite and bounded by to 
Ko — 1, which contradicts the fact that A ^ GDI . Hence Kt < oo a.s. This implies 
that Xf < 1, for all t > 0. Indeed if Xt = 1, for some i > 0, by applying de Finetti's 
Theorem, we deduce that r)t{i) = l,Vz > 1, which contradicts the fact that Kt < oo. 
It remains to show that Xt > for all t > 0. 
For any m > 1, i > 0, we define the event 

= {The m first individuals of type b at time are dead at time t} 

We have 

PiA'P) = (1 - e-'**)"*, 

and then 

p(n„^j") = yt>o. 

From this, we deduce that 3i > 1 such that r]t{i) = 1. The same argument for the 
proof of < 1 shows that Xt > 0, for all t>0. □ 



4.3 The law of X. 



00 



Let X be the proportion of type b individuals at time 0, where < x < 1. As the 
individual at level 1 cannot be pushed to an upper level, we have 

{r/o(l) = 0} C {Xoo = 0}, hence P{X^ = 0) > 1 - x. 
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If Q = 0, is a bounded martingale, so 



P(Xoo = 1) = E(Xoo) = E(Xo) 



If Q > 0, by using (3.5) together with ( |4.2[ ), we deduce that 

P(Xoo = 1) = EXoo < X. 

In this subsection we want to describe those cases where can we decide whether 
P(Xoo = 1) > or P(Xoo = 1) = 0. We first prove 

Proposition 4.4. //AG CDI, then 

P(Xoo = 1) > 0. 

Proof : Since A E CDI, if all individuals at time would be of type b, there would 
be a level J such that the individual sitting on level J at time reaches +oo in finite 
time. Now P(Xoo = 1) > follows from the fact that P(-K"o > J) > 0, where Kq 
denotes the lowest level occupied by a type B individual at time 0. □ 

In the rest of this subsection, we assume that A ^ CDI, and want to decide whether 
F{Xoo = 1) > or P{Xoo = 1) = 0. We shall see that these two situations are possible. 

We first consider the case of a population model dual to Bolthausen-Sznitman's 
coalescent, i.e A{dp) = dp. It is know that in that case the A-coalescent does not come 
down from infinity. 

For any partition with a finite number n > 2 of blocks, the total rate of transitions 
of all kinds in a Bolthausen-Sznitman coalescent equals 

1=2 ^ ^ £=2 ^ ^ 

_ ^ /n\ r(^- l)r(n-^ + l) 

"^vv m (4.5) 

n 

En 

e=2 ^ ' 
= n — 1. 

Recall the definition of Kt. For each n > 1, we have 

P{Ko = n) = (l-x)a;"-\ 

that is Kq follows the geometric distribution with parameter 1 — x. Hence, for each 
X g]0, 1[, E/Cq < oo. We have the following theorem 

Proposition 4.5. If A{dp) = dp, then for any a > 0, 

P{X^ = 1) > 0. 
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Proof : It suffices to show that 



P(Um /Ci = oo) > 0. (4.6) 

t—>-co 

Denote by Ti,T2, . . . the jumps times of the process Kt. Let 

Xji = Kt„ for each n > 1. 

For each integer A; > 1, we define by Vk the size of the jump to the right of the 
individual sitting on level k. Using (4.5), it is not hard to show that conditionally upon 

V \k + Vk, with probability 

I — 1, with probability 

Now, we consider the model where we erase all the arrows pointing to levels above 
Kt for each birth event (t,p) G m (the Poisson point process m is defined by (1.2)). 
We also define Vl the size of the jump to the right of the individual sitting on level 
k in this modified model. We can couple the process X„ with the Markov Chain 1^, 
which jumps at the same birth and death time than Xn, and which evolves as follows 

Y^ = Xo = Ko, 

and conditionally upon Yn = k; 



1 



k + Vlf\ C, with probabihty 
k — 1, with probability 



where C > 2 is a constant to be chosen below. 

Since the "true " model has more arrows than the present model, we have 

Yn <Xn Vn > 1, 

and 

P( lim Yn = oo)< P( lim X„ = oo) = P( lim Kt = oo). 



Using the last statement, to prove ( |4.6[ ), it suffices to show that 

P( lim y„ = oo) > 0. (4.7) 

n— >oo 

For this, we will use Theorem 3.1 in |13]. Let k > 1 denote a fixed integer, and we 
define 

/x(fc) = E(y„+i -Yn\Yn = k) 

^ E(T4^AC)-a 
1 + a 
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and 



= E [{Yn+i - Y^f \Yn = k] 
_ E [{VI, A Cf] + a 



l + a 



Let us compute E(y^ A C). For each p g]0, 1] and k > C, we have 
EM A C) = Ep(Ffc'; Ffc' < C) + CPpiVl > C) 

C k 

= J^(£-l)Pp(Zi + --- + Zfe = £) + C Yl Pp{Zi + --- + Zk = i) 



i=2 
C 



=C+1 



^=2 £=C-\-l 



Next, 



^' Ep(f; a CHdp) = p< - 1) Q ' p<-^(i - + cY: (J) - pf-'dp 



it—n \ / 



k\r{e-i)r{k-£ + i) 



1=2 
C 



m 



k\T{i-l)T{k- l + l) 



m 



k 



i=2 
C 



e=c+i 



£{£ - 1) 



^=1 



We deduce that 



1 ^ 1 



Using the last identity, we obtain 

knjk) _ kE{Vl AC)-ak 
v{k) ~ E{Vl A C)2 + a 
k 



> 
> 



k 



E{Vl AC) -a 



1 

k-1 



{k\ogC - C) - a 
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we choose C such that logC > a, and deduce that 

kfi{k) 
hm — — — = +00. 
v(k) 

From Theorem 3.1 in [13j, the result follows. □ 

In the rest of this subsection, we will give a condition under which P(Xoo = 0). 
Recall the Poisson point measure m defined in the introduction. Let {t,p) be a point 
of the measure m. Let us first give one possible description of the size of the jump to 
the right of an individual sitting on level n at time t~ . 

Let be a random variable with the binomial law B{n,p) and Wq = {L'^ — 1)+. 
If Wq = 0, then there is no jump. If Wq > 1, we put Wq balls in a box Bq. Each 
of those Wq balls, independently of the others and of the value of Wq, is copied with 
probability p, and all copies are put in a box Bi. Let Wi denote the number of balls in 
Bi. Again, each of those is copied with probability all copies being put in a box B2, 
etc.... The just described r.v.'s {W^,i > 1) have the same joint law as the following 
r.v.'s : 



wr 



#{It,pn[l,...,n])-l, ifi = 

#{It,pn[n + l,...,n + W^]), ifi = l 

# (hp n h + Ei=o W^ + l,...,n + Ei=o W^]) ' otherwise 



The size of the jump equals YliLo ^j"- compute 



\i=0 

n-1 



E, 

It is plain that 



'P^o = E (fc + 1) ''P'^'^^ - = np-l + {l-pr. 

k=l ^ ^ 



Bp[wn=pEp[wr^i], i>i. 

Indeed, conditionally upon W^^^ = j, Wj'' is B{j,p) distributed. Consequently 

(^E = (1 - pr'^pm = (1 - p)"' [np - 1 + (1 - pT] . 

Now the mean speed of the movement to the right, starting from n, equals the 
above quantity integrated over u, hence equals 

$(n) = Bp ^E ^" j ''(dp) = _^ V - P)~'[np - 1 + (1 - pTHdp), 
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where v{dp) = p '^A{dp). Using the same argument as in the proof of Lemma 
can easily show that 



4.1 



one 



^ ^ t / ^A(dp) as n t oo. (4.8) 



n Jo p(l-p) 
RecaU that A ^ GDI. We introduce the notation 

The second main result in this section is 

Theorem 4.6. If fx < a, then 

P(Xoo = 1) = E(Xoo) = 0. 

Proof : Suppose fi < a. It follows from lemma 25 of [15] that A ^ GDI. Using 
Theorem 4.3, we deduce that < < 1, for all t > 0. Let Ti be the first time when 
the level 1 is occupied by a type B individual. 

The Theorem is a consequence of the following lemma 

Lemma 4.7. If fj, < a, then 

Ti < (X) a.s 

Proof : Let 

No = inf{n > 1, ^>(n) < a(n - 1)}. 
From (|4.8[) and /x < a, it is plain that Nq < oo. Recall the infinitesimal generator of 



{Kt,t > 0} defined in section 2. Applying (2.1) with the particular choice g{n) = n, 
the process {Mt)t>o given by 



Mt = Kt-Ko- I [^Ks) - a{K, - l)]ds 
Jo 



is a martingale. Let us show that Kt comes back below level No after any time t. 
Suppose that Kq > No and define 

So = mf{t,Kt<No}. 

The process {KtASo)t>o is a positive supermartingale. Consequently, 

KtASo G {iVo,oo} a.s. 

and 

E(i^M5o) < E(ifo) < oo, yt> 0, 

hence by Fatou's lemma, EKqo < oo, consequently K^o = -^o a.s. As the process 
{Kt^So)t>o takes values in the set {A'o, . . . ,oo}, it is easy to deduce that Sq < oo a.s. 
We now define recursively Sk for /c > 1 by the formula 

Sk = mi{t > Sk^i + l:Kt<No} 
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We have 

Pj^^{Ki = 1) > e-<^(^°)(l - e-")^«-^ (4.9) 

Indeed, for {Kt)t>o to reach the level 1 during the time interval [0, 1], it suffices that 
during that time interval there are no arrows between the level 1 and A'^o and during 



the same time all the individuals between 1 and A^o — 1 die. Using (4.9), we deduce by 
the Markov property 



< (1 _ e-<^(^o)(l - e-")^°-^)P(/^Sfe_i > 1) 

< /'i_e-'^(^o)(i_e-")iVo-i^*= 



The result follows by taking the limit as — )• oo. □ 

5 Kingman and A-coalescent 

In this last section we suppose that the measure A is general (i.e A({0}) > 0) and we 
show that the proportion Xt of type b individuals at time t in the population of infinite 
size is a solution of the stochastic differential equation with selection 

Xt = x-a [ Xs{l-Xs)ds 
Jo 

+ I ^/cXs{l-Xs)dB, (5.1) 
Jo 

+ p{lu<x - - Xs-)M{ds,du,dp), 

J[0,t]x]0,l[x]0,l[ 

where c = A({0}),M is the compensated measure M defined in section 3.2, and B 
is a standard Brownian motion. Let {W{ds,du)} be a white noise on (0,oo) x (0,1] 



based on the Lebesgue measure dsdu. We remark that if Xt satisfies (5.1), then Xt is 
a solution in law of the following stochastic differential equation 

Xt = x-a [ Xs{l-Xs)ds 
Jo 

[lu<x^-Xs)W{ds,du) 



+ V~c [ 

J[0,t]x]0,l 

+ / p{lu<x . - Xg-)M{ds,du,dp). 

^fO.tlxlO.lP 



'[0,t]x]0,l 

We first define the model. Recall the process {rjt{i),i > l,t > 0} defined in the 
introduction. The evolution of the population is the same as that described in the case 
A({0}) = except that we superimpose single births, which are described as follows 
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For any 1 < i < j, arrows are placed from i to j according to a rate A({0}) Poisson 
process, independently of the other pairs i' < j' . Suppose there is an arrow from 
i to j at time t. Then a descendent (of the same type) of the individual sitting on 
level i at time t~ occupies the level j at time t, while for any k > j, the individual 
occupying the level k at time t~ is shifted to level k + 1 at time t. In other words, 
Vt{k) = m-{k) for k < j, r]t{j) = Vt{k) = r]t-{k - 1) for k > j. 

By coupling our model with the simplest lookdown model with selection defined in 
[3], it is not hard to show that for N large enough, the individual sitting on level 2N 
at time never visits a level bellow N, that is the evolution within the box (t, i) G 
[0, oo) X {1,2,..., N} is not altered by removing all crosses above 2A''. The process 
^ > l,i > 0} is well-defined. 
For each > 1 and t > 0, denote by X^^ the proportion of type b individuals at 
time t among the first N individuals, i.e. 

1 ^ 

i=l 

Combining the arguments in [3j and section 2.3 (see above), it is easy to show if 
(ryo(i))i>i are exchangeable random variables, then for all t > 0, {rjt{i))i>i are ex- 
changeable. An application of de Finetti's theorem, yields that 

Xt = lim X^ exist a.s. (5-3) 



Using the definition of the model, it is easy to see that 



= X'^ + /Cf + / {K- ' ^' w)M{ds, du, dp, dv, dw) 



Jf0.tlxf0.112 - - 



^ y[0,t]x[0,l]2 

where K,^ is a martingale of jump size ±;^- We have 
Lemma 5.1. 

where, {s) = A(0)Xf (1 - ATf ). 







Proof : For each 1 < i < A^, let be a Poisson process with intensity A(0)(A^ — i). 
At time t £ P^, we have 



N' 



AXi 



N 



if ry^_ (i) = 1 and r/^- (A^) = 
^, if%-(i) = 0andr7i-(Ar) = l 
otherwise. 
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Now, let 



Ai = {i^t{i) = lMN) = ^], 
Bi = {r,t{i) = QMN) = l}. 



We have 

N 

1 

from which, we deduce that 



V{A, I ) = V{B, I ) = -^^X^{1 - XD, 



= A(0)Xf(l-X, 

The result is proved. 
Now, let 

Yf = X^+ [ i}^{X^.,u,p,v,w)M{ds,du,dp,dv,dw) 

J[0,t]x]0,l]4 
J[0,t]x[0,l]2 ~ - 

We have 

Xf = /Cf + y^^, Vi > (5.4) 
From lemma 5.1 , we have VT > 

sup sup I (f^ (s) \< C a.s. 

0<t<TN>l 

Using the last identity, we deduce by Aldous' tightness criterion (see Aldous [1]) that 

{/cf ,t > 0,iV > 1} is tight in D{[0,oo)). 
Since IC^ is tight, there exists a subsequence of the sequence JC^ such that 

/C^ ^ /C weakly in L>([0,oo)), 
where /C is a continuous martingale (since the jumps of JC^ are of size ±;^) such that 

<JC>t= [ cXs{l-Xs)ds, (5.5) 
Jo 

where c = A(0). The main result of this section is 
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Theorem 5.2. Suppose that Xq — )• x a.s, as N ^ oo. Then the [0, 1]— valued process 
{Xt,i ^ 0} defined by (5.3) is the (unique in law) solution of the stochastic differential 
equation 

Xt = x-a [ Xs{l-Xs)ds 
Jo 

[ y/AiO)Xs{l-Xs)dBs 
Jo 

/ p{lu<x ^ - Xg-)M{ds,du,dp), 

J\0,t]x]OAP 



+ 
+ 



(5.6) 



where M is the compensated measure M defined in section 3.2, and B is a standard 
Brownian motion. 

The identification of tlie limiting equation is done similarly as in the proof of The- 
orem 3.6 Strong uniqueness of the solution to (5.6) follows again from Dawson and Li 



[S], and weak uniqueness could also be proved by a duality argument. 

Since Kingman's coalescent comes down from infinity, we have fixation in our new 
model in finite time as soon as A(0) > 0. 
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